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' Invariants of nonlinear gauge transformations of a family of nonlinear Schrodinger 

equations proposed by Doebner and Goldin are used to characterize the behaviour of 
exact solutions of these equations. 

+-> 

3 ; 1 Introduction 

cr. 

In this paper we shall exhibit solutions for members of the family of nonlinear Schrodinger 
equations derived by Doebner and Goldin [jl], |2|. Originally these equations were derived in a 
quantum mechanical context and proposed as possible nonlinear quantum evolution equations. 
There are various contributions in this volume dealing with certain quantum mechanical as- 
pects of this family. E.g. Goldin |J reviews a derivation of the equation and discusses gauge 
transformations of the family, Hennig 0] gives a geometric derivation, Liicke || discusses 
quantum mechanical observables in nonlinear evolution equations and, Malomed || et.al. dis- 
cuss the stability of plain wave solutions, and Mizrahi and Dodonov [^] consider an equation 
generalizing the original equation in |IJ and Gaussian solutions thereof. 

Owing to the presence of many parameters the above mentioned family is very large and 
certain members may be of interest in other fields of physics as well, e.g. PH-||14||. Consequently 
we shall not confine our interest to only certain parameter values (e.g. "small nonlinearity" ) 
but attempt to construct solutions for all values. 

The family of nonlinear Schrodinger equations to be investigated is parameterized by 
the classification parameter D of the unitarily inequivalent group representations and five real 
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parameters D'ci, . . . , D'c^. 

ihd^ = (-^-A + V(x^j ^+i^^ + hD'(f: cjRA^ V ■ (1) 

Here D' also has the dimensions of a diffusion coefficient (so that the Cj are dimensionless) , 
and the nonlinear functionals Rj are complex homogeneous of degree zero, defined by: 

R&] := — , ftW := ^ , R 3 [ii] := ^ , 

*.[*]== *M==^, 

where p := t/^ and J := Im (tp'Vip) = (m/h)j. Our analysis is facilitated if we rewrite the 
family (|T|) wholly in terms of densities and currents. Using the expansion of the Laplacian 
Aip = (iRi[if)] + (l/2)R 2 [ip] - R 3 [i)] - (l/4)R 5 [ip])ip, equation (0) is of the general form 

2 5 

S(u,M) ■= idt*f> ~ i VjRiWpl* ~ Yl - VoVip = (2) 

i=i 3=1 

where u — (z/ 1; u 2 ) and p = (p , . . . , p 5 ) are real parameters and in particular v\ ^ 0. 
Writing if) in the form 

iP = exp(0! + z# 2 ) (3) 
and 9 = (9i, 9 2 ) one finds the following equations: 

AV M (9)i d t 9 l -2u 2 A9 l -u 1 A9 2 -4u 2 (V9 l ) 2 - 2v 1 V9 1 • VB 2 = (4) 

AV M (9) 2 := d t 9 2 +2p 2 A9 x + p,A9 2 + 4(/i 2 + /i 5 )(V^) 2 (5) 

+2(/i! + /i 4 ) V^i • V# 2 + p 3 (V9 2 ) 2 + p V = 0. 

We shall refer to (^) as the amplitude and to as the phase equation. Due to the ambiguity 
of the phase function 9 2 in @ the nonlinear Schrodinger equation (Q) and the amplitude and 
phase equation are not fully equivalent. However, any solution of the amplitude and phase 
equation yields via a solution of the nonlinear Schrodinger equation (0), i.e. AV( y ^){9) = 
0=>S (V)A4) (^) = 0. 

2 Nonlinear gauge transformations 

Before we calculate solutions of the set (|]) of nonlinear partial differential equations (PDEs) 
(fj), let us look for invertible linear transformations of the real functions 0i, 9 2 leaving the 
structure of the amplitude and phase equations (f|,|5|) invariant, i.e. we are looking for a matrix 
A G GL(2), such that 

9>m .= { d 'f^ ) .= ( ^m+AM^t) \ = 
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solve some PDE of the set with parameters (z/,//), if 9 was a solution of (|3],|5|). A short 
calculation shows that the matrix has to be restricted to 

where A ^ 0, hence to a representation of the one dimensional affine group Aff (1). The change 
of parameters according to this transformation is given by 

(z/,/i)^(z/',/i , )=A(A, 7 ).(z/, / i), (8) 

where 

/ v\ , 7 , 

9 

7 7 7 A^3 

fl X = ~JV! // 2 = — z/i - 7^2 - -//i + A/i 2 , ^3 = 7^' ( 9 ) 

2 

7 7 7 

The corresponding transformation of the wave function ■?/> is the one used to linearize a certain 
subfamily of (|I|) in [IS] and partly in [|nj : 

N A{A „0) = V ^ 1+A +^ i( 1 " A +^) . (10) 

Since these transformations leave the probability density p^(x,t) = tp(x,t)if}(x,t) invariant, 
they were called nonlinear gauge transformations ]l7]| (see also || in this volume). Starting 
with a solution of the amplitude and phase equation we can use transformations A G Aff {I) 
to construct solutions of the nonlinear Schrodinger equations of the type (|j) by exploiting the 
implications represented in the following diagram: 

AV M (9) = S> AV AM {A0) = 

^ i/j = exp(9x + i6 2 ) (11) 

S M (?P) = S AM (N A ty)) = 

The change of parameters (||) amounts to the action of the affine group on the eight-parameter 
space, so it is worthwhile to look for a parameterization of (fj) that consists of six parameters 
invariant under the action of Aff (1) and two group parameters. Six functionally independent 
invariants have been proposed in p 



tx = V x fl2 - v-ifi\ , H = fix - 2u 2 , l 3 = 1+ fl 3 / Vx, 1*4 = fa- fiifiz/ v\ , 

(12) 

t5 = Vl{fl2 + 2// 5 ) - Z/ 2 (/ii + 2/i 4 ) + 2v\^j Ux, l = VlflO , 
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and as group parameters we might choose v\ 7^ and Re-expressing the other parameters 
in terms of the new invariants we get 

^2 = 5(^1 - ^) , A*2 = ^r X (2ii - + //?) , 

A*3 = (43 - , A*4 = 44 - A*i + , (13) 

^5 = ^r 1 ( 4 s - ti + t 4 (/xi - 6 2 ) + 2(^1 - ^X^ - 1)) , /X = ^r 1 ^ • 

For the gauge class containing the linear Schrodinger equation ([I]) (D = D' = 0) the 
invariants read 

1 ft 2 

^0 = -7^— , t a = - — ^ , t i = j = 2, ...,5, (14) 
2m 8m 2 

such that linear Schrodinger equations are always represented by invariants of the type to < 

0, t! > 0, 4 2 = • • • = 45 = 0. 

For the calculation of solutions of @ we might limit ourselves to a particular gauge, i.e. to 
a particular choice of the group parameters v%, An appropriate choice is v\ = 1 and /ii = 0, 
and it is sufficient to consider the amplitude and phase equation in this particular "gauge" 

dflx = -12A0! + A9 2 - 2i 2 (V9 1 ) 2 + ■ V<9 2 , (15) 

8 t 9 2 = -2iiA0i -(2i 5 + 24!- 24 4 4 2 - 4^(43 - 1)) (V^i) 2 (16) 
-24 4 V(9! • Vfl 2 - (43 - 1)(V# 2 ) 2 - lqV. 

We will now exhibit various solutions for these equations along with the fact that the "gauge 
invariant" parameters may be used to characterize their different behaviour. 



3 Stationary solutions 

Stationarity is defined as usual by d t p = which is equivalent to d t 9\ = 0. Hence, for stationary 
solutions the amplitude equation ([IB]) reduces to 

- e 2 ) + v0i • vMi - e 2 ) = 0. (17) 

For this equation we have two obvious solutions leading to the following cases (and subsequent 
sub-cases). 



3.1 Plane waves 

Plane wave solutions are obtained if V = in which case setting 9\ (x, t) = const in ([H]) yields 
9 2 (x,t) = k ■ x + k(t) from the amplitude equation such that the resulting phase equation ( |TED 
gives a dispersion relation: 

d t 9 2 {x,t) = (1 - L 3 )k 2 = (1 + ^£^l)e =: u (k) (18) 

n 
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and the solutions thus found are 

ip(x, t) = ipQ exp (i(k ■ x — uj(k))j , ip = const. (19) 

3.2 Nontrivial stationary solutions 

Neglecting constant phase factors the Ansatz 

6 2 {x,t) =i 2 6i{x) -ut (20) 
also solves ( |TTD and leads to the following phase equation 

2n A6 1 +2(n + 6 5 ) (W1) 2 + l V = u . (21) 

3.2.1 Stationary solutions from linear Schrodinger equations 

If 61 7^ 7^ i\ + t 5 and #1 is a solution of (|21|) then the function 



satisfies 



:= exp ( L -l±^ei(x) ) (22) 



2t 2 

" 1 -A<^ + = -</?■ (23) 



From the solutions <p of this equation one obtains the following stationary solutions of the 
nonlinear equation (p]) 

if>(x, t) = (p(x))^ exp (i 1x12 lnta(x)] 2 - iut] . (24) 

V L\ + l 5 J 

As an example we calculate the ground state of the harmonic oscillator in one dimension 

V(x) = fx 2 , for to(ti + is) < 



^')=^(-ife« ,+ 'Hfe«'- W^J) • (25) 

3.2.2 Other stationary solutions 

If L\ = 7^ ig the phase equation reduces to 

2i 5 (Wi) 2 + ^ = w, (26) 

which obviously can only be solved for potentials bounded from above (l > 0) or below (to < 0). 
For such bounded potential separable in Cartesian coordinates V(x, t) = Y^=\ Vj( x j,t) we get 
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or in particular for the free case (V = 0) 9\{x,i) = k ■ x + const and to = 2L$k 2 , hence the 
solution 

ip(x, t) = tpo exp ((1 + ii<2)k ■ x — i2i 5 k 2 tj , ip® = const. (27) 

which is not square integrable. 

If on the other hand L\ + £5 = 7^ L\ the phase equation reduces to 

2l x AOt = to - V , (28) 

i.e. a Poisson equation for 9%, which may be solved by a Poisson integral, if V(x,t) is suitably 
decreasing at infinity 

9i(x,t) =r] X 2 + k-x + C- {2i l )' 1 / G n {x,y)V{y,t)(Ty, 00 = A^n , 

where G n is Green's function for the Poisson equation in n dimensions. Thus in the free case 
the solution 9i(x,t) = k ■ x + V + const leads to the wave function 

ip(x, t) = ipo exp (r/x 2 + k ■ x + 2%ii{r\x 2 + k ■ x) — AiLintj , = const (29) 

which is square integrable only if 77 < 0. 

Finally, if i\ = = 15 then the free phase equation is satisfied by arbitrary 61 with u = 
and we have the following stationary solution 

i/j(x,t) = ipoexp ((1 + 2ii 2 )9i(x)) , ^0 — const. (30) 



4 Gaussian non-stationary solutions 

4.1 Gaussian wave Ansatz 

In the following we shall consider either the case where a harmonic potential V(x) = |x 
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is present or the free case which is thus described by k = [T8|, [Tj|, for a special case pi| . 
Since these cases are separable we restrict our considerations to one space dimension (n = 1). 
Making a Gaussian wave Ansatz for ip is equivalent to making the Ansatze 

(x - s(t)) 2 1 

*iOm) = - l 2a(t y j +^Kt)) (31) 

2 (a:,t) = ^Wa: 2 + 5(t)x + C(t)) . (32) 

Inserting these into the amplitude and phase equation ([15|) and equating equal powers of a; 
yields coupled nonlinear ordinary differential equations for A, B, C, s and a as functions of t. 
This set of equations can be reduced to two ordinary second order equations in a and s 

a = L 3 ^ + 4(i 2 L 3 + L 4 )— + 8 yi „ 5J +2kl (t 33 
a o A u A 

s = Lt + i ^ + ^ s + ^B (34) 

and A, B, and C may be expressed in terms of those. First one has to solve the decoupled 
equation (|33|) and then insert the result into (|34]) and to integrate that equation. 
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4.2 General Gaussian solutions 

For the solution of equations (^) and ( j34[) we will treat the free particle and the harmonic 
oscillator separately. 

In case of the free particle (k = 0) equation (^) can be reduced by the Ansatz 



2 (t) = h(\n(a 2 (t))) . (35) 



dt a 



to the following first order differential equation for h(z): 

2h(z)h'(z) - (1 + i 3 )h(z) 2 - 8(l 2 l 3 + i 4 )h(z) - 32(ii + t s ) = (36) 

Integrating this equation yields in general only an implicit equation for h(z). In certain special 
cases we get an explicit expression, which can be substituted into the Ansatz fl3"5|). As an 
example let us treat the case L\ + 1$ — 0, where 

h(z) = I Cl 6XP ~ 8 ^ f ° r L3 * _1 (37) 

| 4(i 4 - l 2 )z + Ci for l 3 = -1 , 

and we get the implicit solutions 

f*(t) 2xdx 



CT W 2xdx 



'cro 8(i 4 - i 2 ) ln(x) + Ci 
Explicit solutions are obtained for l 3 = 1 



£ — to for i 3 7^ — 1 
t — t for i 3 = — 1 



(38) 



a(f) = ^(2^-^2) + 2(^ + ^4) ^ (39) 

or for i 4 = —l 2 l 3 , t3^1 

cr(t) = (7 ((Tit + l) 1 ^ (40) 

where cri,a 2 are integration constants to be chosen suitably 

For the harmonic oscillator (k > 0) we have not found an analogous Ansatz to reduce fl3"3"|). 
Instead, for t 3 ^ 1 we transform fl3~3|) to an equation for the newly defined function 

q(t) := ((T(t)) 1 - 13 (41) 

such that the following two differential equations for q, s describe the behaviour of the Ansatz 
© in © 

2 ;a±3 
q = A(i 2 L 3 + L^qs-^q - 8(i x + 65) (i 3 - l)g l a-i - 2klq{l 3 - l)q, (42) 

( 2 i 3 g\ 
2(t 2 t 3 + i 4 )g t 3- 1 \s + 2kl s. (43) 

^3 - 1 qj 
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The search for solutions of these equations is facilitated if we view them as one-dimensional 
Newtonian equations with friction. In this way we have the potentials 



U q (q) = 

U s (s) = 
and the friction forces 



-32(ii + t 5 )g - Amoq 2 
-16(ii + 65) lng — 2nioq 2 



( tl +6 5 )( t3 -l) 2 2^±1 
1.3+1 " 



+ Kt (i 3 - l)q 2 



if 13 = 
if L3 = 
if *± 



-3 
-1 

-3, 



-3,-1,1 



KLqS 



F q {q,q) = 4(t 2 t 3 + <, 4 )g'3-ig 

F S (S) = ( 2(t 2 t 3 + i4)g 7 ^ T 



t3 9 



s . 



(44) 
(45) 

(46) 
(47) 



- 1 9, 

Let us treat the case of the harmonic oscillator (/c > 0) first. The asymptotic motion of the 
g-particle depends on the potential U q and the friction force F q . If to < < i\ + t 5 , the potential 
has a minimum. Thus if furthermore 12L3 + ^4 < the friction force F q breaks the motion, and 
q tends to the minimum of the potential and so does a: 

1/4 



lim a(t) = a c 



4(ti + t 5 ) 



lim &(t) = . 



t—*oo 



(48) 



As a consequence at some instant the friction force F s on the s particle becomes negative and 

lim s(t) = , lim s(t) = . (49) 

t^oo t— »oo 

The wave function ip(t) turns out to converge asymptotically to the ground state of the har- 
monic oscillator fl25|). Explicit solutions of (0) with this asymptotic behaviour are quasi-classical 
states with constant width a{t) = cr^ (see section |5~Tp . 

If the above listed conditions are not fulfilled, then the width of the Gaussian wave packet 
might converge to infinity or to zero, i.e. to a constant or a 5 distribution. Furthermore this 
behaviour can depend on the initial condition in the case ti + ts < < to when the potential U q 
has a maximum and quasi-classical states with a(t) = may again be taken to illustrate the 
motion of the wave packet. While the ground state for to < is stable under small perturbations, 
it is unstable for t > 0; a small perturbation leads to a de- or increasing width and a motion 
of the centre s(t) away from the origin. 

For another discussion of the stability of solutions of ([!]) see || in this volume. 



5 Solitary wave solutions 
5.1 Gaussian solitary waves 

Gaussian solitary waves are Gaussian waves whose width remains constant: a(t) = a = const. 
Requiring this to be the case reduces ( |33D to 

K t ( ( To) 4 + 4(t 1 + t5) = 0. (50) 
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In the free case (k = 0) this can only be satisfied if 11 + 15 = 0, for the harmonic oscillator 
cr has to be the width fl4"B|) of the ground state. Equation fl34"|) reduces to an ordinary- 
differential equation of second order: 

s = 2(l 2 l 3 + Li)(yQ 2 s + 2klqs (51) 

For k > and lq < (/o > 0) this is the equation of motion of a damped harmonic (anti- 
oscillator, for k = the equation of motion of a damped free particle. Again the gauge 
invariant expression i 2 i 3 + L 4 determines whether we have damping, pumping or no friction at 
all. 



5.2 Other solitary waves 



Suppose now 9 



(9i(x),9 2 (x,t)) is a solution of the free (V = 0) equations (p0|) and (|21|) 
as in section ^2 and that in addition we have no friction forces for the Gaussian solutions, 
ht>3 + ^4 = 0, then it turns out that 



9i(x, t) = 9\ (x — vt) , 

1 v 2 

9 2 (x,t) = 9 2 (x - vt,t) - -v ■ x + — (1 - L 3 )t 



(52) 
(53) 



is a non-stationary solution of the amplitude and phase equation ftT3[ ) and ftfflj). It should be 
emphasized that the original equation (0) is only Galilei invariant if 13 = = 64 but that the 
transformations (|52|, [53|) take solutions to solutions in the more general case i 2 i 3 + 1 4 = 0. The 
resulting solution of (Q) is then 

ip(x, t) = ipo exp I $i(x — vt) + i(i 2 9 1 (x — vt) — —v ■ x + — (1 — t 3 )t — tot) I , (54) 

where 9\ is a solution of ( pip and ip = const. Consequently the probability density p(x, t) = 
ip(x,t)ip(x,t) becomes a solitary wave 



p(x,t) = |^o 1 2 exp (20i (f - vt)) 

moving with constant speed without changing shape. 

In particular for + t 5 ) < we get square integrable solutions of the type 



ip(x,t) = ipo cosh Ik ■ (x — vt)) tl+ ' 5 exp < i 



2l 2 l 1 
ti + ^5 



In 



cosh (k ■ (x — vt) 



-2- 



ti + ^5 



-JbY 



(55) 



(56) 
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5.3 Solitary waves with arbitrary initial data 

If the parameters are as in the previous section, i.e. t 2 t 3 + t 4 = and in addition ti = = l 5 
then (as was shown in section |3.2.2|) the phase equation (|2lD is satisfied by arbitrary functions 
6i(x). Hence, it follows that for the three parameter family selected by the conditions 

ti = ^3 + U = ^5 = (57) 
the nonlinear equation (0) has solutions 

ip(x,t) = exp (e^x-vt) + i(l 2 9i(x - vt) - \v ■ x + ^-(1 + — )t - ut) ) (58) 

with arbitrary Q\[x). Just as in the previous section this leads to solitary wave motion for the 
"probability density" p(x, t) but now with arbitrary initial data. 



6 Conclusion 

We have presented numerous explicit solutions for members of the family of nonlinear equations 
(fj). These solutions are constructed for the parameters defined by v\ = 1, fi\ — and contain 
the parameter combinations which are invariant under the nonlinear transformation fllCf). Using 
this transformation and the diagram in (|ll]) solutions for all other values of v\ and fii are easily 
obtained. 

The solutions given here contain plane waves, nontrivial stationary solutions, stationary 
solutions associated with a linear Schrodinger equation, Gaussian wave packets, and various 
solitary waves including a three parameter subfamily of (|2|) which permits solitary waves for 
arbitrary initial data. Recently R. Zhdanov has found other explicit solutions of (§) for the 
case 63 7^ ^ l 5 , i\ = t 2 = l± = . 

It seems quite remarkable that despite its nonlinear character so many explicit solutions of 
equation (0) can be found. Apart from its possible use in other fields than quantum mechanics 
the knowledge of ever more explicit solutions may help to shed light on the physical meaning 
(if any) of the various parameters. 
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